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ABSTRACT 



Context. The Crab nebula was observed with the H.E.S.S. stereoscopic Cherenkov-telescope array between October 2003 and January 2005 
for a total of 22.9 hours (after data quality selection). This period of time partly overlapped with the commissioning phase of the experiment; 
Or observations were made with three operational telescopes in late 2003 and with the complete 4 telescope array in January - February 2004 and 
q . October 2004 - January 2005. 

!— I Aims. Observations of the Crab nebula are discussed and used as an example to detail the flux and spectral analysis procedures of H.E.S.S.. 
, The results are used to evaluate the systematic uncertainties in H.E.S.S. flux measurements. 

• Methods. The Crab nebula data are analysed using standard H.E.S.S. analysis procedures, which are described in detail. The flux and spectrum 
^ of y-rays from the source are calculated on run-by-run and monthly time-scales, and a correction is applied for long-term variations in the 
i-H , detector sensitivity. Comparisons of the measured flux and spectrum over the observation period, along with the results from a number of 
>^ ■ different analysis procedures are used to estimate systematic uncertainties in the measurements. 

Results. The data, taken at a range of zenith angles between 45° and 65°, show a clear signal with over 7500 excess events. The energy spectrum 
is found to follow a power law with an exponential cutoff, with photon index T = 2.39 ± 0.03 slal and cutoff energy E c = (14.3 ± 2.1 slal )TeV 
between 440 GeV and 40 TeV. The observed integral flux above 1 TeV is (2.26 ± 0.08 sta t) X 10~"crrr 2 s~'. The estimated systematic error on 
the flux measurement is estimated to be 20%, while the estimated systematic error on the spectral slope is 0.1. 
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1. Introduction 

The Crab supernova remnant (SNR) is an exceptionally well 
studied object, with extensive observations of the system ex- 
isting across the entire accessible spectrum. At a distance of 
2000 parsecs, with an age of 950 years, i t is a prototyp ical 
centre -filled SNR, or plerion, as defined bv lWeiler & Pan agia 
Jl980l) . Within the supernova remnant lies the Crab pulsar, with 
a rotational period of 33 ms and a spin-down luminosity of 
L = 5 x 10 38 erg s _1 . This energy source powers a surrounding 



synchrotron nebula, and polarization me asurem ents exist from 
radio to hard X-ray wavelengths JWilsonlll972l) . indicating the 
non-thermal origin of the radiation detected. The total energy 
available from the pulsar to power the system is of the order of 
10 49 ergs. This is believed to be the power source for produc- 
tion of very high energy (VHE) y-rays. 

The rotational energy of the pulsar is thought to be mostly 
carried away by a relativistic wind of electrons and positrons. 
Interaction of this wind with the surrou nding m edium causes 
a standing termination shock wave jRees & Gunnl 1 19741 
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iKennel & Coronitilll984l). Electron acceleratio n may be due to 
a Fermi-type process ( Achterber g et all 1200 II) or to driven re- 
conne ction of th e alternating magnet ic field at this termination 
shock JCoronitil fl990t iMichell 1 1 994 . The interaction of accel- 
erated electrons with ambient photon fields (in this case mostly 
synchrotron photons) can produce VHE y-rays via the inverse 
Compton process. 

Th e Crab nebul a was discovered at VHE energies in 
1989 JWeekes etall [l989) and emissio n has been con- 
firmed by a number of other experiments. (jGoret et all fl993 : 
Bail lon et all Il993t tM asterson et all 1 19991 lAharonian et al. , 
2000HSmifh et alll2000HAfkins et alil2003l) Due the high flux 
from the source relative to other known TeV sources, and its 
expected flux stability, it is conventionally used as a standard 
reference source for VHE astronomy. At the latitude of the 
H.E.S.S. experiment the Crab nebula culminates at 45°, so ob- 
servations of this source must always be made at large zenith 
angles; the greater effective optical depth of the atmosphere in- 
creases the energy threshold and affects the sensitivity of the 
detector. As the size of the Crab nebula is small compared to 
the H.E.S.S. point-spread function, it may be treated as a point 
source for this analysis. 

A detailed description of the H.E.S.S. detector is given 
here, along with a discussion of the principal sources of sys- 
tematic error in the atmospheric-Cherenkov technique. The 
H.E.S.S. observations of the Crab nebula are then detailed, and 
used as an example in a discussion of the data calibration and 
analysis methods used in source reconstruction and flux and 
spectral measurements. The stability of the y-ray flux and en- 
ergy spectrum is measured using a number of analysis meth- 
ods, and a correction for variations in the long-term optical ef- 
ficiency of the detector is described. Using these results the 
systematic uncertainties on flux measurements with H.E.S.S. 
are quantified. The sensitivity of the detector for source anal- 
ysis is also discussed. It should be noted that while the analy- 
sis presented here is generally used to analyse target s observed 
by H.E .S.S., other techniques are also used, e.g. lLemoine et alJ 

feoodi . 

2. The H.E.S.S. Experiment 

H.E.S.S. is situated in the Khomas highlands of Namibia (23° 
16' 18" South, 16° 30'00" East), at an elevation of 1800 me- 
tres above sea level. The four H.E.S.S. telescopes are placed 
in a square formation with a side length of 120 metres. This 
distance was optimised for maximum sensitivity at the planned 
energy threshold of 100 GeV. 

2.1. The detectors 

The H.E.S.S. telescopes are of steel construction, with alti- 
tude/azimuth mounts capable of precisely tracking any source 
from 0.0 ° to 89.9° i n elevation, with a slew rate of 100° per 
minute (Bolz, 2004J). T he dishes have a Davie s-Cotton style 
hexagonal arrangement JPavies & Cottorl[l957l) with a flat- to- 
flat diameter of 13 m, composed of 382 round mirrors, each 
60 cm in diameter. Thus the effective mirror surface area is 
107 m 2 . Further details of the optical structure are given by 



iBernlohr et alJd2003l) . The mirrors are remotely adjustable un- 
der computer control, and a fully automated procedure is used, 
in conjunction with a CCD camera mounted in each dish, for 
optimal alignment onto the focal plane of each telescope cam- 
era, which is 15 m distant. Due to the rigidity of the dishes, 
this alignment is stable over time scales in excess of one year. 
The stability has been verified by regular monitoring of the 
optical point spread function. Details of the mirror alignment 
system and the opt ical point spread function are discussed by 
ICornilsetalJ(l2003h . 

The H.E.S.S. cameras each consist of a hexagonal array 
of 960 Photonis XP2960 photo -multiplier tubes (PMTs). Each 
tube corresponds to an area of 0.16° in diameter on the sky, 
and is equipped with Winston cones to capture the light which 
would fall in between the PMTs, and also to limit the field of 
view of each PMT in order to minimise background light. The 
camera is of modular de sign, with the PMTs grouped in 60 
drawers of 16 tubes each dVincent etalll2003h . which contain 
the trigger and readout electronics for the tubes, as well as the 
high voltage (HV) supply, control and monitoring electronics. 
The total field of view of the detector is 5° in diameter. 

The trigger system of the H.E.S.S. array consists of three 
levels. First, a single pixel trigger threshold is required, equiva- 
lent to 4 photo-electrons (p.e.) at the PMT cathode within an in- 
terval of 1 .5 nanoseconds. Second, a coincidence of 3 triggered 
pixels is required within a sector - a square group of 64 pixels 
- in order to trigger a camera. Each camera has 38 overlapping 
sectors. Third, when the detector is operating in stereo mode, 
a coincidence of two telescopes triggering within a window of 
(normally) 80 nanoseconds is required. Only cameras which 
have individually triggered are read out in a stereo event. The 
stereo trigger syste m and the trigger b ehavior of the H.E.S.S. 
array is described bv lFunk et al.l (l2004). 

During the first and second level trigger formation, the indi- 
vidual signals from each pixel are stored in two analogue ring 
sampler (ARS) circuits. A high gain and a low gain circuit are 
used to give optimal signal dynamic range. The signals cap- 
tured by each tube are digitised in the drawer, before being col- 
lected by a central CPU in the camera and sent to the central 
data acquisition system (DAQ) by optical ethernet connection 
dBorgmeieretalll2003h. 

The H.E.S.S. experiment commenced observations in June 
2002 with the first telescope. The second telescope was in- 
stalled in February 2003, the hardware level stereo coincidence 
trigger was added in July 2003. In the interim the two tele- 
scopes observed in parallel and an off-line coincidence was 
used for stereo analysis. The third telescope was installed in 
September. The full array was completed in December of 2003 
and has been operational since then. 

2.2. Systematic uncertainties 

The imaging atmospheric-Cherenkov technique depends on a 
form of electromagnetic calorimetry to estimate the energy of 
observed particles. In order to accurately measure the energy 
of the primary particle which gives rise to an air shower, it is 
necessary to understand the relationship between the particle 
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energy and the signal recorded in the cameras. Monte Carlo 
simulations of air showers in the atmosphere are used to predict 
the light yield in the detector, and thus the recorded signal, as a 
function of energy and shower position relative to the observer. 

There are three main contributions to uncertainties in the 
measured air shower information, and thus the absolute flux 
calibration of the detector: 

1. The camera response. The single photo-electron response 
of each PMT varies strongly with the detector voltage and 
is measured using an LED system mounted in front of 
the camera. A second LED system mounted on each dish 
provides a uniform illumination across the camera and is 
used to correct for relative quantum efficiency variations 
of the PMTs, as well as in the reflectivity of the Winston 
cones in front of each PMT. The com plete calibration 
of the H .E.S.S. telescope is described bv lAharonian et all 
ll2004ah . 

2. The optical response of the instrument, including the mir- 
rors, Winston cones, shadowing by structural components 
of the system and the quantum efficiency of the photo- 
cathodes of the PMTs. This response can be measured by 
studying the Cherenkov light from single muons passing 
close to the telescope, assuming the camera response is well 
measured. The use of muons in monitoring t he telescope 
efficie ncies is detailed by iBolzl (2004) and Lerov et al. 
(2004). The optical response of the instrument degrades 
over a timescale of years as, for example, the mirror reflec- 
tivity decreases. This decrease in optical response, relative 
to that used in the Monte Carlo simulations, is taken into 
account in estimating the flux from a y-ray source. This is 
described in detail in section l3~4l 

3. The interactions of particles and light in the atmosphere. 
The atmosphere is the largest and least well understood 
component of a Cherenkov detector, being subject to varia- 
tions in pressure, temperature and humidity. Two important 
effects of variability in the atmosphere are density profile 
variations, which affect directly the height of the shower 
maximum in the atmosphere and thus the intensity of the 
light seen at the telescope, and absorption of Cherenkov 
light in the atmosphere by clouds and dust, which leads di- 
rectly to a reduction in the telescope trigger rates, and in- 
correct y-ray ener gy recons t ructio n. These effects are dis- 
cussed further by Bernlohr (2000). Atmospheric monitor- 
ing devices are used to understand the local conditions un- 
der which the d ata has been recor ded; such measurements 
are discussed bv lAve et alJ l|2£)05). Variations in the atmo- 
sphere can lead to rapid variability in the detector response, 
on a timescale of hours, so runs taken under variable con- 
ditions are rejected, as discussed in section l3~2l 

3. Observations 

3.1. The Crab data set 

The observations of the Crab nebula discussed in this paper are 
summarised in Tabled Data set I was taken between October 
2003 and January 2004 with three telescopes, and 5.3 hours 



of data are included in this sample. These data range in zenith 
angles from 45° to 55°. A further 6.3 hours (data set II) are in- 
cluded from data taken in January 2004 with the complete array 
of four telescopes, at a range of zenith angles from 45° to 65°. 
Data set III includes 11.3 hours, taken between October 2004 
and January of 2005, also with 4 telescopes. All observations 
were taken using a method (wobble mode) whereby the source 
is alternately offset by a small distance within the field of view, 
alternating between 28 minute runs in the p ositive and negativ e 
declination (or right ascension) directions jFominetall ll994). 
This observation mode allows the other side of the field of view, 
which does not contain the source, to be used as a control re- 
gion for estimation of the background level. The wobble offsets 
were varied from 0.5° to 1.5° for the data discussed here. The 
range of time over which the observations have been taken al- 
lows us to study the long-term stability of the H.E.S.S. system. 

Observations of the Crab w ere also taken in 2002 with a 
single telescope, as discussed by Mast erson et al.N2003ll . These 
are not included in this analysis, as the sensitivity is much lower 
than that in stereo mode, and the systematic uncertainties of 
single-telescope observations are greater than for stereo mode. 

3.2. Data quality selection 

Systematic effects on the measured flux and energy spectrum 
(as discussed in section 12. 2> may be ameliorated by rejecting 
data recorded in non-optimal conditions. Although the Crab 
was observed for a total of nearly 45 hours, only 22.9 hours 
are included in this analysis. The remaining observations have 
been rejected as not meeting the run quality criteria. Some were 
made when the sky conditions were less than optimal; the pres- 
ence of clouds or excessive dust in the atmosphere can lead 
to the absorption of Cherenkov light and thus fluctuations in 
the system trigger efficiency, causing systematic uncertainties 
in the measured y-ray flux. Figure^shows the trigger rate as a 
function of time for two runs, one with a stable system trigger 
rate close to the predicted level for this zenith angle, the other 
exhibiting variability due to the presence of clouds. 

Runs for which the mean trigg er rate is less than 70% of 
the predicted value (as discussed bv lFunk et al 1 d2004 l. or for 
which the rms variation in the trigger rate is above 10%, are re- 
jected. As the Crab is mainly visible during the rainy season in 
Namibia when weather conditions are not optimal, the rejection 
rate for these data is not representative of all H.E.S.S. observa- 
tions. The mean system rate is 240 Hz for the four telescope 
data and 180 Hz for the 3 telescope data. 

Quality checks are also routinely carried out in order to re- 
ject runs in which the array tracking system may not be func- 
tioning correctly, leading to errors in the reconstructed posi- 
tion of the source, and thus the flux. The nominal performance 
of the H.E.S.S. tracking system is discussed bv lHofmann et alJ 
d2003l) . The tracking errors reported from the DAQ are moni- 
tored and runs with rms deviations of more than 10 arcseconds 
in altitude or azimuth are rejected, in order to exclude runs in 
which the tracking system malfunctioned. However, no obser- 
vations on the Crab in this study (and in general very few runs) 
have failed this test. As an independent check of the tracking 
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Table 1. Details of the observations of the Crab nebula with H.E.S.S. between October 2003 and January 2005. For the purposes 
of this study, the data has been divided into 3 subsets. 
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Fig. 1. The system rate vs time within a run for a) a 4 telescope run passing the run selection, b) a 3 telescope run failing the run 
selection. The run is rejected due to the instability in the rate caused by clouds passing through the field of view. The difference 
in absolute rate between the two runs is due to the differing zenith angles of the observations, as well as the number of telescopes 
active in each run. 



quality the DC PMT currents JAharonian et al.li2004al) are used 
to estimate the amount of light impingent on each pixel during 
every run as a function of time. A map is then made of the 
sky brightness in the field of view of each telescope. The posi- 
tions of known stars are then correlated with this map, giving 
a measure of the pointing position of each telescope, indepen- 
dent of the tracking system and standard pointing corrections. 
A sample image for the region surrounding the Crab nebula is 
shown in Figure |2 Runs are rejected if the pointing deviation 
is greater than 0. 1°. Again no runs are found to have failed this 
test, which acts as an auxiliary check of the tracking system. 

The presence of bright stars in the field of view may trigger 
the over-current protection of individual PMTs, causing them 
to be turned off. Over-current due to bright star images is pre- 
dicted and the relevant PMTs are turned off in advance, for the 
duration of the star transit through the PMT. Occasionally other 
bright, transient light sources pass through the field of view of 
the telescopes. These are normally meteorites, however light- 
ning, airplanes and satellites may also cause problems. The 
events cannot be predicted and thus trigger the over-current 
protection, causing the PMT to be turned off for the remain- 
der of the run. The mean number of PMTs inactive during each 
run for this and other hardware related reasons is monitored, 
and telescopes with more than 10% of the PMTs missing at 
any one time are rejected from the analysis. The PMT qual- 
ity cut only rejects occasional runs from each data set, and the 



mean number of pixels per telescope disabled during the Crab 
nebula observations is 66, with an rms of 7. 

The observation time after run selection is corrected by tak- 
ing into account the dead-time of the system, when the trigger 
is not sensitive to air showers. Estimation of the live-time (the 
time during an observation in which the telescope system is 
sensiti ve to triggers from the sky) is discussed by iFunk et all 
(12004 . The total live-time for each data set is listed in Tabled 
The systematic error on this estimation is less than 1%. 



3.3. Photoelectron to pixel amplitude calibration 

All of the sources of uncertainty discussed in section l2~2l af- 
fect the pixel intensities, and thus the energy reconstruction of 
each shower. The trigger efficiency and energy threshold are 
also affected, contributing to uncertainty in absolute flux mea- 
surements. Part of this effect, namely the camera response, is 
continuously calibrated using the LED system; the accuracy of 
this calibration system can be evaluated by comparing the re- 
constructed pixel intensities from Monte Carlo simulations of 
air showers with similar measurements on real data. Figure [3] 
compares the intensities in the pixel with the third highest sig- 
nal amplitude, as seen in real data for all triggering events, with 
those from Monte Carlo simulations of protons (which com- 
prise a large majority of events detected) at the same zenith 
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Fig. 2. Sky map of DC pixel current (in units of ji/A, measured 
in a window of 16 yus) for the region surrounding the Crab neb- 
ula. The positi ons of stars (read from the Hipparcos catalogue 
- |Perrvrnan&RSANl997ln match well with the peaks in the 
map, which correspond to peaks in sky brightness. The B and 
V magnitudes of each star are also given for comparison. The 
brightest star, £ Tauri (B magnitude 2.8), has caused the camera 
to disable PMTs, thus there is no current measure for the region 
close to this star. 



angle. As the trigger requirement is 3 pixels above a certain 
amplitude, the amplitude distribution of the 3rd pixel shows the 
effective pixel trigger level of each telescope. The peak in the 
distribution is determined by the trigger threshold of the tele- 
scope and the energy spectrum of the simulated and real events. 
It can be seen that the Monte Carlo simulations reproduce well 
the data Ot 2 /dof= 15.5/23). 

3.4. Optical efficiency correction 

The optical efficiencies of the H.E.S.S. telescopes change over 
time, probable causes include degradation in the mirror reflec- 
tivity. This effect happens on a timescale of years and can be 
monitored using images of local muons in the field of view, for 
which the light yield can be predicted, as discussed in section 
l2~2l 

The effect manifests itself as a reduction in the image in- 
tensity for each event, compared to the intensity expected from 
Monte Carlo simulations. This causes a shift in the absolute en- 
ergy scale of the detector, as events are reconstructed with ener- 
gies which are too low. This effect is corrected by incorporating 
a scaling factor into the energy estimation for each event. The 
image intensity used in the energy estimation is scaled by the 
ratio of the mean optical efficiency (over the telescopes) for the 
run (Eff r „„) to the mean optical efficiency as derived from the 
Monte Carlo simulations (Eff, m . ). The corrected energy is then 
derived from this scaled image intensity in the standard man- 

Eff 

ner. The distribution of the relative optical efficiency ( gff "" ) f° r 




10' 
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Fig. 3. Distribution of the intensity of the 3rd highest pixel for 
4 telescopes (points), compared with Monte Carlo simulations 
(solid lines). The overall flux of the simulated data has been 
adjusted to fit the real data. 
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Fig. 4. Distribution of the relative optical efficiency per run for 
data sets I (solid line), II (dashed line) and III (dotted line). It 
can be seen that the relative optical efficiency is significantly 
decreased in data set III relative to the other data sets. 



the runs included in this analysis is shown in Figure^ The ap- 
plication of the optical efficiency correction in flux estimation 
is discussed further in section ROl 
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Fig. 5. Definition of simple Hillas parameters, calculated for 
a y-ray image, which may be approximated as an ellipse. 
Important parameters for this analysis are the width, length, 
distance. An image from a second telescope is superimposed 
to demonstrate the geometrical technique for source position 
reconstruction. The parameter 8, which is the magnitude of the 
angular offset in shower direction reconstruction, is also shown. 
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Fig. 6. Distribution of excess events in 9 2 for the complete Crab 
data set, after event selection and background subtraction. The 
Monte Carlo derived point-spread function described in equa- 
tion 1 is also shown, normalised to the excess distribution. The 
vertical lines denote the 6 2 selection cuts listed in Table|2 



4.2. Stereo reconstruction 



4. Analysis 

After a set of images of an air shower has been recorded, they 
are processed to measure Hi llas paramet ers based on the sec- 
ond moments of the image ilHillasl F985). These parameters 
are then used for event selection and reconstruction. A diagram 
illustrating the parameter definitions is shown in Figure|5] 



4.1. Image cleaning and moment analysis 

The first step in the moment analysis procedure is image clean- 
ing. This is required in order to select only the pixels contain- 
ing Cherenkov light in an image. Other pixels, which contain 
mainly night sky background (NSB) light are not used in the 
analysis. Images are cleaned using a two-level filter, requiring 
pixels in the image to be above a lower threshold of 5 p.e. and 
to have a neighbour above 10 p.e., and vice versa. Cleaning 
thresholds of 4 p.e. and 7 p.e. have also been shown to work 
satisfactorily, but may be more sensitive to uncertainties due 
to NSB light variations. This method selects spatially corre- 
lated features in the image, which correspond to air shower 
Cherenkov light. This method tends to smooth out shower fluc- 
tuations in a simple and repeatable manner. 

After image cleaning, an image of a y-ray shower approx- 
imates a narrow elliptical shape, while images of background 
hadronic events are wider and more uneven. The Hillas param- 
eters are then calculated for each cleaned image; these parame- 
ters are the basis for event selection. The total amplitude of the 
image after cleaning is also calculated, along with the mean 
position of the image in the camera, which corresponds to the 
centroid of the ellipse. 



The arrival direction of each event is reconstructed by trac- 
ing the projected direction of the shower in the field of view 
(which corresponds to the major axis of the image) to the point 
of origin of the particle. For stereo observations it is possible 
to intersect the major axes of the shower images in multiple 
cameras, providing a simple geometric method of accurately 
measuring the shower direction; more details, including meth- 
ods t o further imp r ove th e reconstruction accuracy are given 
by [Hofmann et al. (1999), method I from that paper is used 
here. Images are only used in the stereo reconstruction if they 
pass the selection cuts on distance (to avoid camera-edge ef- 
fects) and image intensity. If less than two telescope images 
pass these cuts the event is rejected. 

Figure [6] shows the excess distribution of 6 2 for data sets 
I-III, including events with two, three and four telescopes; 6 
is defined in Figure [5] it is the angular offset between the re- 
constructed shower direction and the true direction of the Crab 
nebula. The distribution of reconstructed shower directions is 
usually expressed in units of 9 2 , as this ensures a constant solid 
angle on the sky per bin. The value of the cut on reconstructed 
shower direction is thus given in units of degrees 2 in Table|2]for 
various sets of selection cuts, and plotted in Figure[6] A strong 
excess is seen close to zero, corresponding to events coming 
from the direction of the Crab nebula. This distribution defines 
the accuracy in the reconstructed arrival directions for y-ray 
events from a point source and is described by the point spread 
function (PSF). This function can be approximated by the sum 
of two, one-dimensional Gaussian functions: 



PSF = A 



exp 



2cr 2 ) 



+ A r 



•i exp 



Q2\ 



2a- 2 , 



(1) 
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This function is fitted to the 9 2 distribution for simulated 
Monte Carlo y-rays. For simulations at 50° zenith angle the 
standard deviation parameters cr 1 and cr 2 are 0.046° and 0.12° 
respectively. The relative amplitude of the second Gaussian, 
A re / is 0.15, while the absolute amplitude (A) is proportional 
to the number of events in the fit. The 68% containment radius 
for 50° zenith angle is 0.12°, while that for 10° zenith angle is 
0.10°. The point-spread function is shown in Figure|6]with the 
amplitude parameter (A) fit to the Crab nebula data. The^- 2 /dof 
of this fit is 53/29. 

The position of the centre of the Cherenkov light pool, 
which corresponds to the projected impact point of the origi- 
nal particle track on the ground, can also be reconstructed by 
intersecting shower axes, projected into the plane perpendicu- 
lar to the system observing direction. It is vital to reconstruct 
this position in order to accurately measure the amount of light 
originally emitted by the shower and thus the shower energy. 
The rms error on the reconstructed impact parameter, which 
is the projected distance of the extrapolated shower track to a 
telescope, for Monte Carlo simulations is less than 10 metres 
for events falling within 200 metres of the centre of the array. 

4.3. Scaled parameter analysis 

The mean scaled widt h method, similar to that used by the 
HEGRA collaboration JDaum et all fl997). is used to classify 
images as either y-ray like or hadron like, in order to reject non 
y-ray background events. The main difference to the HEGRA 
method is in the definition of the scaled parameter itself, in the 
HEGRA case this is defined as p TC . = p/(p). A lookup table 
is used to predict the mean width and length for a y-ray as a 
function of the amplitude of the shower image in the camera 
and impact parameter. Then the value for a particular event (p) 
can be compared with the expected value (p) according to the 
formula: 

Psc = (P - <P»/cr p (2) 

The mean value (p) and the scatter cr p for an event vary 
with the image amplitude and impact distance, as well as the 
zenith angle. Lookup tables are generated for 13 zenith angles 
(Z) from 0° to 70°, based on Monte Carlo simulations. The true 
impact parameter of the simulated shower is used in filling the 
table. 

When analysing real data, the reconstructed impact param- 
eter is used along with the image amplitude for each telescope 
image to find (p) and <x p in the lookup table. Linear interpo- 
lation (in cos(Z)) between the two nearest simulated values is 
then done to find the correct value for a particular observation 
zenith angle. The mean reduced scaled width (MRSW) and the 
mean reduced scaled length (MRSL) are then calculated by av- 
eraging over the telescope images passing the image amplitude 
selection cut for each event: MRSW = (2 te i Psc) /N te i. 

Figure0(a) shows a comparison between the MRSW from 
Monte Carlo simulations of protons and y-rays and from real 
data at a zenith angle of 50°. It can be seen that the data (be- 
fore selection cuts) correspond well to Monte Carlo simulated 
protons, as expected, while there is good separation between 



the data and Monte Carlo simulated y-rays, which are chosen 
to have a photon index (T) of 2.59, similar to that p reviously 
measured for the Crab nebula iAharonian et all EoOO). 

4.4. Selection cuts 

Selection cuts on the mean scaled parameters, image intensity 
and 6 2 are simultaneously optim ized to maximi se the detec- 
tion significance (<x, as defined bv lLi & Mai <1983lfl for sources 
with typical fluxes and energy spectra. The optimisation pop- 
ulation consists of a mixture of y-ray simulations (selected to 
give the desired flux and spectrum for optimisation) and real 
background data. In the presence of background, the signifi- 
cance achieved for a given source increases with the square 
root of the observation time; instrument performance is there- 
fore characterised by crj -^f (hr _I ). The optimised cuts yield the 

maximum crj J?(hr _1 ) for a source of that type. It should be 
noted that the optimum selection cuts in any analysis depend 
on the energy spectrum of the Monte Carlo simulations used in 
the optimisation procedure, and it may be necessary to optimise 
selection criteria separately for much harder or much softer 
energy spectra. As a rule this is not done however in source 
searches, in order to preserve the a priori nature of the analy- 
sis. The background data used in the optimisation is then not 
further used, in order to avoid the possibility of optimising on 
background fluctuations and compromising the statistical in- 
dependence of the results. The selection cuts are summarized 
in Table |2] A number of alternative sets of cuts are presented 
here, which are used for analysis of H.E.S.S. sources. It will be 
shown that the reconstructed flux and spectrum of the Crab are 
consistent for these various selection criteria. 

The standard set of selection cuts has been optimised to 

give the maximum <x/ Jt (hr -1 ) for a flux 10% of the Crab 
(standard cuts), with a similar spectrum. The Hard cuts are 
optimised for a source with a flux 1 % of the Crab flux, and a 
r of 2.0. These cuts give a higher significance for weak, hard 
spectrum sources, at the expense of energy threshold and cut 
efficiency. The hard cuts are also useful as they reduce the sys- 
tematic uncertainties in sky-map reconstruction by reducing the 
numbers of background events, relative to the signal. They also 
give a narrower PSF than the standard cuts, as the higher inten- 
sity cut selects better reconstructed events. A set of Loose cuts 
have been also optimised to give the maximum significance for 
a strong source, similar to the Crab, and a T of 3.0. The lower 
intensity cut here reduces the energy threshold of the analysis 
relative to the standard cut, and the fraction events passing the 
cuts is higher. When conducting source searches, the standard 
cuts are always used unless there is an a priori reason to expect 
a very hard or very soft spectrum from the source. 

For analysis of large extended sources the cut on 6 2 is usu- 
ally set to be larger than the extension of the source, so that 
effectively all y-rays from the source pass this cut. In order to 
demonstrate the effect of this strategy, a version of the standard 
cuts is described, with the 2 cut set to a much larger value. 
These are referred to as extended cuts for the purposes of this 
paper. It should be emphasised that only the standard selection 
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Fig. 7. The distributions of mean reduced scaled width (MRSW) for Monte Carlo y-ray simulations (F = 2.59) a) with Monte 
Carlo proton simulations (r = 2.70) and actual off data before selection cuts, b) with real events from the direction of the Crab 
nebula (data set II) after selection cuts and background subtraction. All distributions are for zenith angle = 50°. The vertical lines 
denote the standard cuts described in Table|2] 



cuts are used in searches for point sources, extended source 
searches are carried out using an a priori 6 2 cut suited to the 
source size, and trials are taken into account when testing mul- 
tiple source extensions. 

Figure shows the distributions of MRSW, after standard 
selection cuts and background subtraction (see following sec- 
tion), for real data and simulations (with the same mean zenith 
angle). The cut on MRSW is not applied for this plot. The stan- 
dard MRSW selection cuts are indicated, it can be seen that 
the cuts select y-ray-like events. The small shift between the 
data and Monte Carlo simulations seen in this plot is due to 
differences in the optical efficiency; simulations with reduced 
efficiency (as in Figure 0} agree well with the data. It can be 
seen that this shift has a negligible effect on the efficiency of 
the scaled parameter cuts. 

5. Signal extraction and background estimation 

When estimating the flux of y-rays from a particular direc- 
tion in the sky it is necessary to estimate the background level, 
due to non y-ray events with directions reconstructed close to 
the source direction. The significance of the excess after back- 
ground subtraction is the n determined usi ng the likelihood ra- 
tio method described by Li & Ma| dl983h . For the purpose of 
background estimation the distribution of background events is 
usually assumed to be azimuthally symmetric within the cam- 
era field of view. However, zenith angle dependent effects or 
variations in the NSB level across the field of view may in- 
troduce non-radial variations in the background level. A radial 
profile of the relative rate of background events passing shape 
cuts (the background acceptance) in the field of view is shown 
in Figure|8](dashed line). For comparison, a number of test ob- 



servations (duration 30 minutes each) were made at a range of 
offsets from 0° to 2.5° on the Crab nebula. It can be seen that 
the relative rate of excess y-ray events passing cuts (points) fol- 
low the background acceptance closely out to 1.5° offset in the 
camera. The Monte Carlo predicted y-ray rate for this zenith 
angle is also shown (solid line), this is described in section l6"3l 
When estimating the background, first the reconstructed 
shower direction for each y-like event (i.e. an event that passed 
the shape cuts) is filled in a two dimensional histogram (so- 
called raw sky-map). The on signal for a given point in the sky 
is determined by selecting events within a circle around that 
point with radius 6 cut . Two techniques are used to derive esti- 
mates of the background level within this region of the field of 
view, and are described below. 

5.1. Reflected background model 

The simplest background estimation technique uses the signal 
at a single position in the raw sky-map, offset in the opposite 
direction relative to the centre of the field of view, to estimate 
the background. This technique is used in the standard wobble 
observation mode (described in section^. However, it suffers 
from relatively low statistics in the measurement of the back- 
ground level due to the choice of a single reflected background 
position, as well as possible systematic effects caused by local 
inhomogeneities at the background position. 

The generalised reflected background technique, which is 
also suitable for wobble mode, uses a number of background 
regions equidistant from the observation position, as illustrated 
in Figure|9] The combined events from these positions are used 
to estimate the background at the on position, scaled by the rel- 
ative area of the on and off regions. In the case of a larger on 
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(degrees^) 
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-2.0 


2.0 


-2.0 


0.9 


0.0125 


80 


2.0 


Hard 


-2.0 


2.0 


-2.0 


0.7 


0.01 


200 


2.0 


Loose 


-2.0 


2.0 


-2.0 


1.2 


0.04 


40 


2.0 


Extended 


-2.0 


2.0 


-2.0 


0.9 


0.16 


80 


2.0 



Table 2. Selection cuts optimised for various purposes, as described in the text. Cuts are applied on MRS W and MRSL, as well as 
on the distance (6) from the reconstructed shower position to the source. A minimum of two telescopes passing the per-telescope 
cuts, on image amplitude and distance from the centre of the field of view, are also required. Standard cuts, as well as hard, loose 
and extended cuts, as described in the text, are listed. 
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Fig. 8. A comparison of the predicted relative y-ray rates 
(from MC simulations) and those measured from data (from 
4-telescope test observations on the Crab nebula), as a function 
of off-axis angle in the field of view. Also plotted is the relative 
radial acceptance for background events passing selection cuts. 
This agrees well with the relative y-ray acceptance out to 1.5°. 




-0.5 



X(°) 



Fig. 9. Schematic illustrating the background regions described 
in the text. The observation position of the telescopes is marked 
by a cross, while the target position is marked by an X. The 
on region surrounding the target position is marked by a cross 
hatched circle. The ring-background region is filled by hori- 
zontal lines, while the refiected-background regions, at con- 
stant offset from the observation position, are filled by diagonal 
lines. The two background regions have equal area in this case. 



integration region the number of background regions is reduced 
to eliminate overlapping. The normalisation, a, is the ratio of 
the solid angles of the on and off regions. As the off positions 
are the same distance from the centre of the field of view as the 
on position, correction for the relative radial background ac- 
ceptance of the detector is not required. However, this method 
cannot be used for positions closer to the centre of the field of 
view than the radius of the on region, as the background posi- 
tions would overlap with the source position. As all of the data 
described here is taken in wobble mode, this method is used in 
this analysis for flux and spectral measurements. 

The refiected-background method may also be susceptible 
to systematic effects caused by non-radial variations in the ac- 
ceptance, especially for large offset positions in the field of 
view. Non-radial effects are strongest for data with only two 
telescopes, where the trigger efficiency can vary with the az- 



imuth direction of the shower impact point on the ground rel- 
ative to the telescope positions. This could be corrected by ap- 
plying a small non-radial correction term across the field of 
view. This correction has not been applied in the analysis de- 
scribed here. However, in order to minimise systematic effects 
due to non-radial acceptance variations, the direction of the 
wobble offset is generally alternated run by run on either side 
of the target position. 

5.2. Ring-background model 

The ring-background technique determines the background 
for each position in the field of view using the ba ckground 
rate c ontained in a ring around that position ( Piihl hofer et all 
2003). The internal and external radii of the ring are here cho- 
sen such that the ratio of the areas of the off to on regions is 
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(a) (b) 

Fig. 10. a) Uncorrelated 2-d plot of y-ray excess from the Crab nebula, the reconstructed sky-map in RADec J2000 coordinates 
around the source position is shown, (b) The distribution of significance per bin in an uncorrelated significance sky map for the 
same data. The open circles denote the distribution for points in the map within 0.4° of the Crab position, while the filled circles 
represent points further from the source. The Gaussian function, fitted to the second distribution, has a mean of 0.04 + 0.006 stat 
and a standard deviation of 0.98 ± 0.004 st!it . A mean of 0.0 and standard deviation of 1.0 is expected for an unbiased significance 
distribution of the background. 



close to 7, which makes for a convenient compromise between 
area within the ring and distance from the on position. The 
inner ring radius is chosen to be significantly larger than the 
on region, in order to avoid signal leakage into the off region. 
The normalisation (a) is given by the area ratio modified by a 
weight factor to account for the radial background acceptance 
in the camera. The ring around the on position is illustrated in 
Figure [9] When estimating the background for a test position 
close to a known source like the Crab nebula, the source posi- 
tion is cut out of the background ring in order to avoid signal 
pollution in the off region for the test position. This method 
has the advantage of allowing background estimation for all 
positions in the field of view. However, since the number of 
events at positions surrounding the source is used to estimate 
the background in the direction of the source, it is most suitable 
for sources with an small angular extent relative to the field of 
view of the detector. Figure [TOf a) shows an excess map of the 
sky in the vicinity of the Crab nebula, after background sub- 
traction. Figure [Tor b) shows the distribution of significance of 
the excesses in each bin in the sky map. It can be seen that the 
significance is distributed normally in the off source regions of 
the map (filled circles), while the region close to the Crab neb- 
ula (open circles) shows a significant excess. 

The ring-background method is less suitable for spectral 
analysis of sources than the reflected method as the background 
acceptance may not be constant as a function of energy, thus 
the background level may not be correctly estimated for the 



entire energy range of the spectral analysis. Thus this method 
is not used for the main spectral and flux analysis in this paper, 
however results are derived using this method for comparison 
purposes. 

6. Energy reconstruction and effective areas 

6.1. Energy reconstruction 

The energy of the primary particle of a y-ray shower is esti- 
mated for each telescope as a function of the image amplitude 
and impact parameter using a lookup table. The lookup table 
contains the mean energy for Monte Carlo y-ray simulations 
as a function of total image amplitude and the simulated true 
impact parameter. As for the scaled parameters, the lookup ta- 
bles are created for a number of zenith angles and the resulting 
energy is estimated by linear interpolation in cos(Z), and aver- 
aged over the triggered telescopes for each event. Events with 
relative error in the reconstructed impact parameter greater than 
25% are not used in the lookup table creation, in order to reject 
poorly reconstructed events, which may bias the lookup table. 
Events with a distance greater than 2° from the centre of the 
field of view are also rejected. 

In the case where no estimate is available in the lookup 
table for a particular event, due to a lack of Monte Carlo statis- 
tics at that combination of impact distance and total amplitude, 
an alternative lookup is used with coarser impact distance bin- 
ning. This occurs on average for 0.3% of events. In the case 
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Fig. 11. Relative bias ((E reco - E, nw ) /E, nw ) in estimation of en- 
ergy as a function of energy for three zenith angles. The vertical 
lines represent the safe energy thresholds for spectral analysis 
at each zenith angle. 
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Fig. 12. The distribution of the relative error in the recon- 
structed energy per event for Monte Carlo simulated y-rays 
with a power law energy distribution (above 440 GeV) with 
F = 2.6 at 50° zenith angle. The root mean square (rms) width 
of this distribution is 16%. The width of the fitted Gaussian 
distribution is 14%. 



where the optical correction is applied, the energy is estimated 
using the corrected image intensity, as described in section l3T4l 
The energy estimated for each telescope is averaged to give the 
mean energy for an event: E est = (£ te i E te i) /N te i. 

6.2. Energy resolution 

The error in the reconstructed energy (Ae) for a particular sim- 
ulated y-ray event with true energy E true and reconstructed en- 
ergy E reco is defined as Ae = (E reco - E true )/E true . The mean 
value of Ae is shown as a function of E, rue in Figure^J F° r en- 
ergies close to the threshold, there is a bias due to a selection ef- 
fect, whereby events with energies reconstructed with too high 
a value are selected. In order to make an accurate energy spec- 
trum it is necessary to define the useful energy range, so as to 
avoid the region of large energy bias. First the lowest energy 
bin in the bias histogram with a bias of less than 10% is found. 
The useful lower energy threshold is the maximum energy of 
this bin plus 10%. This safe energy threshold is also indicated 
in Figure \^\ for each zenith angle. This energy threshold for 
the analysis is increased to take account of the shift in the en- 
ergy scale when the optical efficiency correction is applied. It 
can be seen that the energy bias above the safe threshold does 
not depend on the zenith angle, up to energies in excess of 60 
TeV. 

The distribution of Ae, for y-rays simulated with a power 
law spectrum with F = 2.6, at 50° zenith angle, is shown in 
Figure^]for the standard analysis. The energy resolution for 
a particular energy range is defined as the width of this distri- 
bution. Events in this plot are selected above the safe threshold 



of 0.44 TeV in order to avoid the effect of the energy bias. The 
energy resolution defines the optimum binning for spectral re- 
construction, as well as defining the minimum energy width of 
any resolvable spectral structure. It is possible to improve the 
energy resolution slightly by selecting only those events with 
higher telescope multiplicities and smaller impact parameters, 
at the expense of reduced event statistics. 

6.3. Effective areas 

The y-ray flux from a source is estimated from the number of 
excess events passing the selection cuts for a particular data set 
using the effective area of the instrument. The effective area is 
a function of the zenith angle and offset of the source from the 
pointing direction, the energy of the event and the particular 
selection cuts used. The effective area is modeled from Monte 
Carlo simulations by counting the fraction of simulated events 
which trigger the detector and pass the selection cuts. This ef- 
fective area has been estimated in two ways: as a function of 
the Monte Carlo energy of the simulated events (Atrue) and as 
a function of the reconstructed energy (Areco). While ktrue 
does not depend on the energy spectrum of the simulated y- 
rays, the finite energy resolution makes Areco sensitive to this; 
the effective areas are usually estimated assuming a power law 
distribution of F = 2.0. However, when estimating a flux using 
events binned in reconstructed energy, it is correct to use Areco 
to estimate the effective area for each bin, when estimating the 
integrated effective area over the whole energy range one may 
use Atrue. In order to avoid a bias in the spectral reconstruction 
when using Atrue, it is necessary for the energy spectrum of the 
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Fig. 13. The effective collecting area of the full H.E.S.S. array versus energy a) as a function of true Monte Carlo energy and b) 
as a function of reconstructed energy for observations at zenith angles of 20°, 45°, and 65°. The vertical lines denote the safe 
energy thresholds for each zenith angle, increasing in zenith angle from left to right. 



simulations to match that of the data, this is discussed further 
below. The effective areas as a function of true energy and re- 
constructed energy (for the standard selection cuts) are shown 
in Figure^] for three zenith angles. As the effective area of 
the telescope system depends strongly on the zenith angle of 
the observations, it is determined for a range of angles and the 
value for a particular energy and zenith an gle is de termined 
by line ar interpolation in log(£') and cos(Z) (Aharonian et al., 
Il999al) . 

In order to simplify the application of the optical correc- 
tion discussed in section 13.41 this correction is not applied in 
estimation of the effective area for each event. Since the dis- 
tribution of Cherenkov light scales with the shower energy to 
a good approximation, the detection probability and hence the 
effective area depends primarily on the amount of light arriving 
at the camera, and not on the absolute energy of the y-ray. Thus 
it is not necessary to recalculate the effective area lookup tables 
when the actual optical efficiency changes by a small amount. 
For larger changes this correction breaks down due to effects 
of the system trigger on events near the energy threshold. The 
selection cuts are made using the image intensity without op- 
tical correction, and the corrected intensity is only applied in 
the energy estimation. This method has been tested on Monte 
Carlo simulated sources with reduced optical efficiency, and it 
was verified that the correct flux is reconstructed. 

Instead of using the safe energy threshold as introduced 
above, the energy threshold for a set of observations has also 
been common ly defined as the peak i n the differential rate vs. 
energy curve (Kononel ko et all 1 1999). This is formed by fold- 
ing the effective area curve, as plotted in Figure fOf a) with the 



expected y-ray flux from the source. This energy threshold is 
generally slightly lower than the safe threshold defined above, 
which is designed to ensure an accurate spectral reconstruction. 
There may even be a significant y-ray signal below the safe 
threshold. The vertical lines in Figure define the safe en- 
ergy thresholds for each zenith angle. Figure[2]shows the pre- 
dicted peak-rate energy threshold and y-ray rate for a Crab-like 
source, based on simulations and projection f rom the Crab flux 
as m easured by the HEGRA collaboration d Aharonian et all 
2000). Table |3] gives the pre-cut energy threshold as a func- 
tion of zenith angle, along with the equivalent after the various 
selection cuts. It should be noted that the energy threshold, by 
this definition, depends on the spectrum of the source. 



The effective area also varies with the position of the source 
in the field of view of the instrument. As larger energy show- 
ers are preferentially detected at higher impact distances, they 
appear closer to the edge of the field of view. Thus truncation 
of images for sources closer to the edge of the field of view 
tends to reject events at higher energies. Monte Carlo simula- 
tions are made at a range of source offset positions in order to 
make effective area curves and the effective area is interpolated 
for a particular observation. Table|3]also shows the event rates 
before selection cuts and after each of the selection cuts de- 
scribed above, as a function of zenith angle for a source similar 
in flux to the Crab nebula. These predictions are based on the 
effective areas estimated for a source offset by 0.5° from the 
observation position. The cut efficiencies in each case are also 
shown, with the (peak rate) energy thresholds. 
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Table 3. y-ray rate predictions from simulations for the standard, hard and loose selection cuts. This table is valid for a source 
with an energy spectrum similar to the Crab, for observations at an offset of 0.5° (the usual observing mode). The cut selection 
efficiencies and (peak rate) energy thresholds in each case are also given. 
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Fig. 14. The peak-rate energy threshold before and after selec- 
tion cuts versus zenith angle, for three sets of selection cuts as 
described in Table|5] The safe threshold for spectral analysis in 
each case is slightly higher. 



7. Flux and spectral measurements 

The data from observations of the Crab have been analysed us- 
ing the technique described, individually by data set as outlined 
in Tabled and combined. The analysis has been carried out for 
each of the sets of selection cuts as outlined in Table |2] using 
the reflected background method with 5 background regions. 
Table 0] outlines the numbers of events passing cuts in the on 
and off regions, as well as the background normalisation (a), 
the number of excess events, the significance (<r) of the excess, 

the rate of y-rays passing cuts and the crj -yjt (hr _1 ). Data Sets 



I-III are combined to give a total result. For example, the mean 
rate of y-rays for the standard selection cuts is 6.0 y min 1 

with a significance of 27 cr/ ^jt (hr _1 ). It can be seen that the 
rate of events passing cuts is strongly dependent on the zenith 
angle of the observations, as well as on the selection cuts used. 
Comparisons between Tableland Table[3]show that the effec- 
tive area estimation correctly reproduces the y-ray rate for the 
various data sets (within statistical errors). The mean rate for 
data set III, at a mean zenith angle of 54°, is predicted to be 
4.8 y min -1 , while the measured value is (4.9 + 0.1 stat )y min -1 . 
The loose cuts at 55° zenith angle keep 53% of the possible y- 
rays, making them more suitable for spectral studies of strong 
sources. 

7.1. Run by run flux measurements 

For the purpose of producing a light curve of the y-ray flux 
from a source, the integrated flux above the threshold en- 
ergy is calculated for each time period {t s , art to t st0 p), assum- 
ing a particular spectral form for the source, such as a power 
law with photon index F and flux normalisation Io in units 
of cm -V^eV -1 . The excess number of events seen from a 
source (6) is given by the following: 

5=\ l (—)- r Atrue(E,Z(t))dtdE (3) 

Here Atrue is the effective area as a function of zenith an- 
gle Z and true energy E from Monte Carlo simulations. The 
flux normalisation can then be calculated by integrating the ef- 
fective area up to some upper cutoff energy E f and over the 
integration time. The value of E c is imposed by the range of 
the Monte Carlo simulations and is normally above 100 TeV. 
The integral flux above threshold is usually quoted above the 
threshold energy for the observations, or alternatively above 1 
TeV. Here the latter convention is used for simplicity. 

The fluxes measured for each data set included in this anal- 
ysis, along with the mean fluxes for data sets I-III are given in 
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Table 4. Events passing cuts in on and off regions for the Crab, listed by data set along with excesses, significance and y-ray rates. 
Various selection cuts described in Tableware also compared for data sets I-III combined. The results using the ring-background 
model are given (denoted as Ring). The integrated flux from the Crab above 1 TeV is shown also, as described in section lTTl The 
X 2 /dof for a fit to a constant flux for the data set is given, as is the percentage run-by-run rms deviation in the flux. 
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Fig. 15. The run-by-run light curve of the integral flux above 1 TeV for a) data set I b) data set II c) data set III. All efficiency 
corrections as discussed in the text have been applied to these data. 



Table 0] The mean flux is also given for the various sets of se- 
lection cuts described in Table[2] The rms variation in the mean 
flux per data set is 8%, while the rms variation in the run-by- 
run fluxes is 14.9%, the typical statistical error on a single run 
is 5% for moderate zenith angles and offsets. Figure IT5l shows 
light-curves for data sets I, II and III, while the distribution of 
run-by-run flux is shown in Figure^] for all of the data. 

It can be seen that the long-term variations in the run-by-run 
flux after the correction for changes in the detector optical effi- 
ciency are small compared to the short-term variations, mainly 



due to atmospheric effects. The mean flux in data sets II and 
III, taken in January 2004, is 10% higher than that seen in data 
set I, which was taken in October of 2003. This difference is 
smaller than the rms spread of either data set, and can be ex- 
plained by differences in the atmospheric conditions between 
the two periods. No correction is made to the run-by-run flux 
for short term variations in the atmospheric conditions. Such 
corrections are under study and will be the focus of a future 
paper. 



F. Aharonian et al.: Observations of the Crab Nebula with H.E.S.S. 



15 




2.5 3 
Flux> 1TeV(x10" 11 ' 



Fig. 16. Distribution of run by run fluxes for data sets I-III. The 
fitted Gaussian distribution has a mean of (2.21 + 0.06 stat ) x 
lO^'curV 1 andao-of(3.58+0.6 stat )x lO^cirrV 1 . 

Systematic errors in the flux due to selection cuts and effec- 
tive area estimates have been studied by applying the various 
selection cuts described to the data and measuring the integral 
flux. The results are outlined in Table @] The rms in the mea- 
sured fluxes is 15%. A possible systematic error in the integral 
flux measurement due to the background estimation method 
has been tested by applying the ring -background method, as 
described above to calculate a flux. It can be seen that the re- 
constructed flux differs only slightly with the two methods. 

7.2. Energy spectrum 

The energy spectrum of the Crab has been measured using 
the data described here. The method used for deriving t he en- 
ergy spectrum is similar to thos e de scribed bv iMohantv et all 
ill 9981) . lAharonian et al.ldl999bl) and lAharonian et al.1 d2004bl 
An energy spectrum is fit for each data set and for the combined 
data. 

The bin size for the energy spectrum is set depending on 
the overall significance of the signal. The maximum possible 
energy bin is defined by the simulations, which extend to above 
400 TeV at 50° zenith angle. Only energies above the safe 
threshold as defined above are used for spectral determination. 

In each energy bin ; above the minimum energy the differ- 
ential flux is calculated by summing over the on source events 
N on , weighted by the inverse of the effective area (Areco) as 
a function of the reconstructed energy of each event. The nor- 
malised sum of the weighted off events N /f is then subtracted. 
The difference is weighted by the live-time for that bin (T) and 
the bin width (AE,): 
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In the case where runs are combined with varying zenith 
angles, and thus varying useful energy thresholds, the live-time 



Fig. 17. The ratio of the reconstructed effective area (estimated 
assuming an photon index of 2.0) to the true effective area per 
energy bin, for true photon indices from 1.1 to 3.2, based on 
Monte Carlo simulations at 45° zenith angle. 



is calculated for each energy bin separately. The result is scaled 
by the appropriate live-time to give the y-ray flux in each en- 
ergy bin. The error on the flux in each energy bin is estimated 
using standard error propagation. A spectral energy function 
(for example a power law distribution) is fit to the flux points 
using the least-squares method. The maximum energy for spec- 
tral fitting is chosen so as to have a significance in that bin 
greater than 2<x. 

Since the estimation of the effective area as a function of 
reconstructed energy, Areco, depends on an assumed spec- 
tral slope, it is strictly correct to adjust Areco for the fit- 
ted spectrum and then re -fit, repeating until the fit converges 
( Ahar onian etafl. fT999b). Figure [TTl shows the bias as a func- 
tion of energy introduced by using an assumed photon index 
of 2.0 for the effective area estimation, given true spectra with 
indices ranging from 1.1 to 3.2. It can be seen that at 440 GeV, 
the differential flux for a source with an intrinsic photon index 
of 2.6 is overestimated by 5%, while the differential flux given 
a true photon index of 1.5 is underestimated by 4%. For ener- 
gies well above threshold the bias is less than 5% for a wide 
range of photon indices. Thus the effect of this correction on 
the Crab spectrum is small and was neglected for this analysis. 

The number of excess events and significance is given for 
each energy bin along with the differential flux in Table[5] Only 
statistical errors are given here. A significant signal is seen in 
every energy bin from the threshold energy of 440 GeV up to 20 
TeV, and a marginal signal is seen at the maximum bin of mean 
energy 30.5 TeV. The spectrum is shown in figure ED The fit 
of a power law function to the combined data with the standard 
analysis cuts yields T = 2.63 + 0.02 stilt and differential flux nor- 
malisation at 1 TeV Io = (3.45 +0.05 stat ) x lO^cirrV'TeV- 1 
(PL in Table|6). 

In the combined spectral fit there is evidence for a steepen- 
ing of the energy spectrum; a fit of a power law with an expo- 
nential cutoff: H = Io (E/lTeV) ( - r) exp(-E/E c ), gives a dif- 
ferential flux normalisation at 1 TeV Iq = (3.76 + 0.07 stat ) x 



lO^'cirrV'TeV 1 , with F = 2.39 + 0.03 stat and a cutoff en- 
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Fig. 18. a) Energy spectra for data sets I (filled circle), II (open circles) and III (filled triangle), fit residuals to the common power 
law fit are also shown. The dashed line indicates the best fit power law spectrum, while the solid line denotes the fit including 
an exponential cutoff, b) Combined average energy spectrum for data sets I-III, fit residuals to the combined power law fit are 
shown. 
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Table 5. Flux measurements for each energy bin in the combined spectral fit on data sets I-III, as plo tted in Figure IT8T b). The 
flux errors are error-propagated 68% Feldman-Cousins confidence intervals jFeldman & Cousinslll998l) . 



ergy E c = (14.3 + 2.1 stat ) TeV. The^ 2 for this fit is 15.9 with 
nine degrees of freedom. This compares with a^ 2 of 104 with 
ten degrees of freedom for the straight power law fit, thus the 
fit including an exponential cutoff is clearly favoured. A broken 
power law fit (BPL in Table[6j: 

f=Mtr( 1+ ( t )" s ) !,ri - r " 



gives a differential flux normalisation Io = (3.43 +0.07 slal ) x 
lO^'cirrV'TeV-', with T = 2.51 + 0.02 slal below the break 
energy of E e = (7.0 + 0.1 slat ) TeV. The photon index above the 
break is F = 3.3 + 1.5 sta t, and the^ 2 /dof of the broken power 
law fit is 28.6/8. This fit includes a term (S) for the width of the 
transition region, which is fixed to 0.3. 
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Table 6. Flux and spectral measurements of the Crab, divided up by data set as outlined in Tabled f° r a power law fit with an 
exponential cutoff. Results for the various selection cuts described in Tableware also compared. The results for a power-law fit 
(PL) and for a broken power-law fit (BPL) are also given. The spectral fit estimated using the ring background model is given 
(ring), as is that estimated without the optical efficiency correction (uncorr.). Only statistical results are show n in the table. S imilar 
measurements fr om other experiments a re given for comparison, the Whipple results is taken from iMohantv et all Jl998h . the 
CAT results from Master son et alJJl999l) and the HEGRA results from Aharonian et al. (20043)- 



The flux and spectral measurements for the separate data 
sets are summarised in Table [6] The rms spread of the photon 
index from data set to data set is 0.04. The rms spread in the 
integral flux, calculated from the fitted spectrum, is 15%, with 
a statistical error of 2% on the integrated flux for the combined 
data. Figure flTjl a) shows the energy spectral points superim- 
posed for the three data sets used in the spectral analysis, the 
residuals about the combined fitted spectrum are shown under- 
neath. Figure [T8l b) shows the energy spectral fit, along with 
residuals, for the combined data sets I-III. The spectral fits for 
the various selection cuts are also included in Table|6| as is the 
fit for the ring background model. For comparison the fit spec- 
trum calculated without the optical efficiency correction is also 
included. 

Figure ^] shows a 2-d plot of the fitted photon index F 
against the flux normalisation for each data set analysed, here 
a simple power law fit has been made for simplicity. The er- 
ror contours, estimated using the least-squares method, are also 
shown. It can be seen that the three data sets are compatible at 
the 2<t level; data set II includes large zenith angle data, and 
is fitted with a slightly softer spectrum, caused by the higher 
energy threshold of these observations and the significant cur- 
vature seen in the spectral measurements. 

The combined data sets I-III have also been analysed with 
the various selection cuts described in Table|2] and a spectrum 
fitted. The rms spread in the photon index between the various 
analyses is 0.08, however this includes the effect of the very 
different energy threshold for the hard selection cuts, which 
may give rise a softer photon index if the source spectrum is 
intrinsically curved. The rms spread of the reconstructed inte- 
gral flux is 8%, which indicates that the reconstructed flux is 
not strongly dependent on the details of the analysis method. 
The use of the ring-background method results in a flux and 




Fig. 19. Contour plot of the x 2 fit error as a function of the 
power law parameters Fo and F for data sets I (filled circle), 
II (open circle) and III (filled triangle, dotted lines). The 68%, 
95% 99.9% error contours are shown, and The best fit spectral 
parameters are marked in each case. 

spectral slope similar to that reconstructed using the standard 
reflection method. 

7.3. Estimation of systematic errors 

The systematic error on the absolute flux is estimated from the 
various independent contributing factors, as discussed in sec- 
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Table 7. Summary table showing the various estimated contri- 
butions to the systematic flux error 



tion l2.2l These errors are summarised in TableQ The total es- 
timated systematic error, after correction for degradation in the 
optical efficiency, is 20%. The sources of systematic error in- 
clude uncertainties due to the shower interaction model and the 
atmospheric model used in the Monte Carlo simulations. Also 
included is the estimated uncertainty in the flux due to the effect 
of missing pixels, which has been conservatively estimated at 
5%, and the effect of uncertainty in the live-time measurement, 
which is less than 1%. The uncertainty in flux due to selection 
cuts is estimated from the rms of the flux and spectral slope 
measurements detailed in Table [6] as is the uncertainty due to 
the background model. The run-by-run rms over the entire set 
of data is 15%, this is thought to mainly be due to variations 
in the atmosphere, and thus is included in the systematic error 
as an independent factor. The rms of the spectral estimations 
for the various datasets in Tablets used to estimate the uncer- 
tainty in the spectral slope, which is 0. 1 . 

8. Conclusions 

A strong signal has been detected from the Crab nebula during 
the commissioning phase of the H.E.S.S. instrument and with 
the complete instrument. An energy spectrum has been mea- 
sured, with a differential spectrum described by a power law 
with slope T = 2.39 + 0.03 st!it + 0.09 sys and an exponential cutoff 
at (14.3 + 2.1 stat + 2.8 sys ) TeV. The integral flux above 1 TeV is 
(2.26 + 0.08 stal + 0.45 sys ) x lO^cnrV 1 . 

Marginal steepening in the spectrum measured on the 
Crab nebula has been previously claimed by Ahar onian et alJ 
(2004b) in studies of the Crab with the HEGRA experiment. 
Figure |20| compares the Crab spectrum from this study with 
measurements by HEGRA, CAT and Whipple. Acceptable 
agreement is seen up to 10 TeV between the experiments, al- 
though the CAT result gives a somewhat steeper spectrum; the 
rms variation in integral flux between the four experiments is 
22%. Above 10 TeV the energy spectrum as seen by H.E.S.S. 
steepens significantly, in particular compared to HEGRA. 

The softening seen in the Crab spectrum at high ener- 
gies is consistent with models of inverse Compton emission 
due to a population of electrons extending up to PeV ener- 
gies. Due to the high magnetic field in the Crab nebula, the 
dominant target photon field for emission is probably created 
by synchrotron em ission from the same electron population 
jHillas et all F998>. More detailed models of y-ray emission in 
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Fig. 20. Comparison between the spectral measurements of 
CAT (open squares), Whipple (open circles) and HEGRA 
(open triangles) and the results if this study (filled circles). The 
residuals are shown for each spectrum relative to a power law 
fit to the H.E.S.S. data (dotted line). 



the Crab nebula are discus sed bv lAtovan & Aharonian! Jl996): 
Ide Jaeer & HardindJl992l) . 

Given the agreement between the Monte Carlo simulations 
and the data, one can use the simulations to predict the time re- 
quired to detect a point source of a certain strength as a function 
of zenith angle. In Figure|^we show the time as a function of 
signal strength required for a 5cr detection at 20° zenith angle, 
for the selection cuts described here. The H.E.S.S. array is ca- 
pable of detecting a point source with a flux of 1 % of the Crab 
nebula in 25 hours, or alternatively detecting a source of similar 
strength to the Crab in 30 seconds. The sensitivity for extended 
sources decreases approximately linearly with the source ex- 
tension. This sensitivity is unprecedented in the field of VHE 
astrophysics and opens a new window for sensitive and precise 
measurements of VHE y-ray sources. 
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Fig. 21. Sensitivity of the HESS array, expressed as the amount 
of time required to detect a signal at the 5<x level, as a function 
of the flux of the source, for a source of similar spectral slope 
to the Crab nebula (for observations at 20° zenith angle, 0.5° 
offset from the source). A minimum of 10 excess events is also 
required for very short observations. Shown are curves for the 
standard, hard and loose selection cuts. 
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